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Abstract 

A critical appraisal is given of a recent analysis of the quark-mass and finite-size dependence of 
unquenched lattice QCD data for the nucleon mass. We use this forum to estimate the boundary 
of the chiral regime for nucleon properties. 



1 



How low in quark masses must lattice QCD simulations of baryon properties go in order 
to reach the regime where chiral perturbation theory (xPT) is an extrapolation tool with 
controlled errors? In interesting recent work 0, 0] unquenched lattice QCD data for the 
nucleon mass at large pion masses, > 500 MeV, has been analyzed using baryon xPT 



to 0{m'^) (without including the A as an explicit degree of freedom). These papers fit the 
quark-mass dependence of unquenched lattice QCD data for the nucleon mass (including the 
physical point) with natural values of the strong-interaction parameters and then successfully 
predict (in Ref. [3]) the finite-size dependence. These results are quite remarkable. However, 
little attempt is made in Refs. jH I4I to gauge the uncertainties associated with the quark- 
mass extrapolation procedure. Given the provocative results that are found and the rather 
large quark masses that are simulated, it is essential to perform an error analysis. We will 
therefore do so in this note. And, in so doing, we will provide an estimate of the quark-mass 
boundary of the chiral regime for nucleon properties jsf . 

The nucleon mass at 0{mVj in the chiral expansion, evaluated using infrared regulariza- 
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Following Refs. we choose the parameter set gA = 1.267, / = 131 MeV, C2 = 3.2 GeV~^ 
and C3 = —3.4 GeV~^. One then easily fits the remaining low-energy constants to un- 
quenched lattice QCD data 0, over a wide range of pion masses ^Ij,!^]. For instance, 
in Table lU (Fit 1 (IR)) we list a set of parameter values ^ which give rise to the solid curve in 
Fig. ^ (left panel). While one sees an excellent description of the data with natural values of 
the strong-interaction parameters (which are consistent with n — N phase shift analyses), one 
must quantify the errors associated with omitted higher-order effects in order to estimate 
the reliability of perturbation theory over the range of fit pion masses. 



TABLE I: Fit parameters. The scale-dependent parameter ei is evaluated at = 1 GeV. 





/ (MeV) 


gA 


Mo (GeV) 


ci (GeV-i) 


ei (GeV-3) 


Fit 1 (IR) 


131 


1.267 


0.880 


-0.95 


2.44 


Fit 2 (IR) 


124 


1.520 


0.872 


-1.22 


4.90 


Fit 1 (DR) 


131 


1.267 


0.885 


-0.92 


3.61 


Fit 2 (DR) 


124 


1.520 


0.879 


-1.07 


6.79 



There are many ways to estimate the reliability of the chiral expansion and Refs. [l|, |2| 
suggest one robust method. The difference between a matrix element at a given order in the 
chiral expansion evaluated with the physical values of qa and / and evaluated with their 
chiral-limit values g^^ and Z^^-* is a measure of the importance of higher-order effects. Ref. [3] 
estimate g^^^ ~ 1.2 and Z^^-* ~ 124 MeV and find little variation in the fit curve. However, 



We fit the full expression in eq. while Refs. |lll2| fit to an expanded form; hence the small differences 
in the values of the parameters. 



2 




FIG. 1: Left panel: The nucleon mass in baryon xPT (computed using infrared regularization) 
at 0{m'^) (parameters given in TahleU (Fit 1 (IR)))vs. the pion mass squared. Symbols are 
unquenched lattice QCD data taken from Ref. S Q/- The solid curve is the fit curve and 

the dashed curves give maximal variations due to the spread in chiral limit values of f and gA, 
as explained in the text. The gray region corresponds to the error associated with omitted higher 
orders. Right panel: Same but with a second fit, (parameters given in Tabled (Fit 2 (IR))) with 
corresponding spread in the chiral limit values of f and gA- (Note that there are two curves going 
through the data.) 



this is no surprise as these chiral-hmit values leave the ratio qa/ f unchanged. Here we will 
consider how variation of the chiral-limit parameters affects the chiral expansion. 

It is worth considering what is known experimentally abouto?'* and f'^^K At one-loop 
order in the chiral expansion one has the well-known formulas js], | 
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where di% and I4 are (scale- independent) low-energy constants ^. 

The parameter I a is determined from K^^ and K12 decays and from a two-loop analysis 
of vr — TT scattering For simplicity, we will follow Ref. j9| and fix /^^^ = 124 MeV, keeping 
in mind that not accounting for variation of f^^^ will necessarily lead to an underestimate 
of the final error. An analysis [1, 0| of the process vrA^ — > vrvrA^ at one-loop order in 
the chiral expansion provides three distinct determinations of di^. —0.91 ± 0.74, —1.01 ± 
0.72 and —1.76 ± 0.85 GeV~^. Taking the Icr limits of the three determinations gives 
-2.61 GeV"^ < rfi6 < -0.17, or = 1.42 ± 0.10. A recent analysis suggests that the A 
(which is not included as an explicit degree of freedom in Refs. 



the determination of die and suggests a lower central value: g\ 



(0) 



plays a crucial role in 
1.20 ±0.10 Q . Taking 



^ The barred constants contain chiral logarithms and therefore the formulas in eq. 151) are useful only at 
the physical value of the pion mass. Extrapolation formulas are available in Refs. [lo[[l]| . 
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—2.61 GeV~^ < die < 2.43 GeV~^ encompasses both analyses, i.e. Ag^^^ = 1.10 — 1.52. As 
naturalness suggests l^iel ~ 1, this range of values does not introduce anomalously large 
low-energy constants into the chiral expansion. This experimental/theoretical uncertainty 
therefore provides an estimate of the importance of neglected higher orders in the chiral 
expansion ^. In Fig. [T] (left panel), we illustrate (dashed curves) the spread of the fit curve 
when one replaces / and qa by their chiral limit values, including the range of rfig values. 

It is important to realize that all curves encompassed by the gray region of Fig. Q (left 
panel), between the dashed curves, differ only by terms that are higher order, (9(mg/^), in 
the chiral expansion. The spread of curves in Fig. [T] (left panel) therefore suggests a 50% 
error associated with neglected higher orders in the chiral expansion at the lattice point 
with the lowest pion mass. We also perform a second fit with select chiral-limit values of / 
and Qa, see Table H] (Fit 2 (IR)), and again considering variation in the chiral limit values ^. 
The first and second fit are shown in Fig. ^ (right panel), together with the spread (dashed 
and dotted curves) of the fit curves when one replaces / and qa by their chiral limit values, 
including the range of rfig values. The spread of curves in Fig. [T] (right panel), suggests an 
80% error at the lattice point with the lowest pion mass. 
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FIG. 2: Left panel: The nucleon mass in haryon xPT (computed using dim reg with MS J at 0{ml) 
(parameters given in Tahle^ (Fit 1 JDR)J) vs. the pion mass squared. Symbols are unquenched 
lattice QCD data taken from Ref. 0,0, Q;^/- The solid curve is the fit curve and the dashed curves 
give maximal variations due to the spread in chiral limit values of f and gA, cls explained in the 
text. The gray region corresponds to the error associated with omitted higher orders. Right panel: 
Same hut with a second fit (parameters given in Tabled (Fit 2 (DR))) with corresponding spread 
in the chiral limit values of f and gA- 



^ We emphasize that even if dig (and therefore ) were known with high precision, varying Jig over a range 
of natural values would continue to be a legitimate way of estimating errors due to omitted higher-order 
effects. 

^ Notice that although the parameter ei(l GeV) appears rather large, we have defined ei as in Refs. 
which has absorbed a factor of 4 into its definition as compared to the parameter to which naturalness 
arguments should be applied 15]. Hence, all values of ei in TableQ]are technically quite natural. 
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It is interesting to compare the error analysis performed using infrared regularization 
with the same analysis performed using dimensional regularization (dim reg) with MS; as 
physics must be independent of the regulator, so must the results of the error analysis ^ 
The nucleon mass in dim reg with MS at 0{mg) in the chiral expansion is given by 



Mn = Mo- Aciml - j^ml 
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< . (3) 



Again one easily finds a fit to unquenched lattice QCD data over a wide range of pion masses. 
In Table m (Fit 1 (DR)) we list a set of parameter values which give rise to the solid curve in 
Fig. 121 (left panel). Again we consider variations in the chiral limit values of / and qa (gray 
region) and perform another fit (see Table |l] (Fit 2 (DR)) and Fig. |21 (right panel)). We find 
very little difference in the error analyses for the two regularization schemes. 

Our results are summarized in Fig. |21 which plots the errors in the nucleon mass extrap- 
olation curves taken from Figs. Q and El the curves take into account the error associated 
with the gray regions in Figs. Q and El Clearly the chiral expansion is nicely convergent 
in the vicinity of the physical pion mass. In our view the dashed curves (right panels of 
Figs. H and EI) give a conservative estimate of the errors. If one is willing to tolerate a 20% 
error associated with neglected higher-order terms for an (9(m^) calculation, then the error 
analysis would indicate that one is in the chiral regime for 300 MeV. If one instead 
takes the errors given by the solid curves (left panels of Figs. [Hand El), willingness to tolerate 
a 20% error would indicate that one is in the chiral regime for m,r^ 400 MeV. It may appear 
odd that our estimate of the boundary of the chiral regime is less than the kaon mass, which 
governs the convergence of SU (3) baryon xPT. We stress that the convergence of the chiral 
expansion is process and flavor dependent; for instance, some observables in SU{3) baryon 
XPT converge well and others do not. 

The method used here to estimate errors is, of course, only one method among many to 
quantify the reliability of perturbation theory. For instance, by comparing various orders in 
the chiral expansion, Ref. finds that one is in the chiral regime for m.^^ 500 MeV. An 
interesting finding of this work is that the 0{m'^) correction remains a small perturbation 
on the (9(mg/^) result for ttItt < 600 MeV. It is important to emphasize that error analyses 
such as that presented in this note and in Ref. [17| are merely indicative. We believe that 
we have given a conservative estimate of the errors, as they presently stand. Naively, one 
way of weakening the strength of higher-dimensional operators, and thereby reducing the 
error, is to include the A as an explicit degree of freedom in xPT ^. It may also be possible 
to reduce the extrapolation error at larger pion masses by computing M^r, and / in the 
same lattice simulation in order to fix the lattice values of (iig^ and Z4. It is clear from the 
results presented here that a definitive determination of the boundary of the chiral regime 



^ Refs. P, Q use the nomenclature "relativistic xPT" to describe infrared regularization; we dislike this 
appellation as it may be misinterpreted to suggest that there is physics in the regulator. 

^ As very- few quantities have been computed including the A at non-trivial orders in baryon xPT, it is 
presently unclear whether or not this is the case. 

^ This parameter also provides a major source of uncertainty in the quark-mass dependence of the deuteron 
binding energy [liJ, lllj 
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FIG. 3: Errors in the nucleon mass extrapolation curves vs. the pion mass squared. The solid lines 
correspond to the errors abstracted from the left panels of Figs. Q] and\^ they each arise from the 
effect of Ag^^ and f^^^ on a single fit (Fit 1), performed with the physical values of gA and f. 
The dashed lines correspond to the errors abstracted from the right panels of Figs.^\and\^- same 
as above but supplemented with the effect o/Ag^^ and f^^^ on a second fit (Fit 2), performed with 
select chiral-limit values of gA and f. The horizontal dotted lines indicate 10% and 20% errors. 
The star indicates the physical pion mass. 



will not be possible until these strong-interaction parameters are well determined. Finally, 
it is clear that by imposing various prejudices on our model-independent analysis, one may 
shrink the gray regions of Figs. Hand|21to any desired size. 

In this note we have provided evidence that the quark-mass dependence — and by associa- 
tion, the finite-size and lattice-spacing dependence — of currently-extant unquenched lattice 
QCD data for the nucleon mass are not presently described by a perturbative effective field 
theory with a controlled error estimate. Hence these lattice data cannot be reliably extrap- 
olated to predict nucleon properties. In our view, a symbiotic relationship between lattice 
QCD and baryon chiral perturbation theory which will lead to first-principles predictions of 
baryon properties with meaningful errors must await smaller lattice quark masses. 
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AC05-84ER40150, under which the Southeastern Universities Research Association (SURA) 
operates the Thomas Jefferson National Accelerator Facility. 
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